Geometric-phase interference in a Mni2 single-molecule magnet with four-fold rotational symmetry 
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We study the magnetic relaxation rate F of the single-molecule magnet Mni2-tBuAc as a function of mag- 
netic field component Ht transverse to the molecule's easy axis. When the spin is near a magnetic quantum 
tunneling resonance, we find that F increases abruptly at certain values of Ht- These increases are observed just 
beyond values of Ht at which a geometric-phase interference effect suppresses tunneling between two excited 
energy levels. The effect is washed out by rotating Ht away from the spin's hard axis, thereby suppressing the 
interference effect. Detailed numerical calculations of F using the known spin Hamiltonian accurately repro- 
duce the observed behavior. These results are the first experimental evidence for geometric-phase interference 
in a single-molecule magnet with true four-fold symmetry. 
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PACS numbers: 75.45.+j, 75.50.Xx, 03.65.Vf 

Geometric -phase effects are responsible for many fascinat- 
ing phenomena in classical and quantum physics from how a 
cat rights itself while in free fall to the dynamics of charged 
particles in electromagnetic fields, e.g. the Aharonov-Bohm 
effect [1]. One formulation of geometric-phase effects in- 
volves a path-integral approach in which the interference of 
paths is modulated by the geometric phase difference between 
the paths |2-5]. Such interference effects can reveal the under- 
lying symmetries of the system's Hamiltonian. The dynamics 
of spins provide a natural way to explore quantum geometric 
phases for, as Berry showed in his pioneering work IQ], a sys- 
tem near a degeneracy point can be mapped onto a spin in an 
effective magnetic field. 

Single-molecule magnets (SMMs) are an ideal test bed for 
exploring spin geometric-phase interference. In these sys- 
tems, each molecule behaves as a large spin with a well- 
defined Hamiltonian determined by the symmetry of the 
molecule and its environment [7]. Although there are many 
molecules in a crystal, the interactions between them are typ- 
ically weak and the sample behaves as an ensemble of nom- 
inally identical large-spin objects. In most SMMs of inter- 
est, the spins have a large anisotropy barrier separating the 
preferred "up" and "down" directions. This leads to hys- 
teresis and slow relaxation between these easy-axis direc- 
tions. A geometric -phase effect can lead to interference be- 
tween tunneling paths, thus modulating the rate at which 
spins flip direction. In a ground-breaking experiment, Werns- 
dorfer and Sessoli |Q] found oscillations in the probabil- 
ity of magnetization tunneling as a field applied along the 
hard axis modulated the interference between two tunneling 
paths. This observation confirmed a theoretical prediction by 
Garg |4] and ignited intense theoretical study of related phe- 
nomena |9- 17]. Geometric -phase interference between tun- 
neling paths has been observed in other SMMs that have ef- 
fective two-fold symmetry, where tunneling involves the inter- 



ference between two equal-amplitude paths IIISll . It has also 
been seen in antiferromagnetic molecular wheels lfl9ll and in 
SMM dimers i2d l2lll . However, such interference effects 
have not been seen in the bellwether SMM Mni2Ac, pre- 
sumably because its nominal four-fold symmetry is broken 
by solvent disorder i22ll . Here we report the observation of 
a geometric -phase interference effect in [Mni20i2(02CCH2- 
*Bu)i6(CH30H)4] CH30H (hereafter Mni2-tBuAc), a vari- 
ant of Mni2Ac that is free of solvent disorder and main- 
tains its four-fold rotational symmetry 123142711 . Unlike pre- 
vious observations of geometric -phase interference, which 
involved ground-state tunneling, the interference effect de- 
scribed herein is observed in the thermally assisted tunneling 
regime where tunneling takes place near the top of the barrier 
The interference effect provides a fingerprint that affords an 
unprecedented ability to clearly identify which levels partici- 
pate in the thermally assisted process. 

In 2002, Park and Garg UM and, independently, Kim Q, 
predicted that an interference effect should be observed in 
SMMs with fourth-order transverse anisotropy, described by 
the Hamiltonian 
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-DS^, + (C/2){Sl + St)-gfiBS-H, (1) 



where S± ~ Sx ^ iSy. In zero field, such a spin sys- 
tem has the classical energy landscape shown in the inset of 
Fig. [T] Applying a magnetic field Ht along one of the four 
hard directions (± x and ± y for C > 0) preserves reflec- 
tion symmetry through the z-hard plane, allowing for interfer- 
ence between equal amplitude tunneling paths that virtually 
pass through the saddle points in the landscape. This inter- 
ference induces oscillations in the tunneling probability as a 
function of Ht- Mni2-tBuAc is reasonably well described by 
the above Hamiltonian with the addition of sixth-order terms 
consistent with four-fold symmetry t25l] : 



2 



n = 



-DSl - ASt - A' si + (C/2) {Si 
+ (C"/2) (5f (5t + ^l)+/..c.) 
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where £> = 0.568 K, A = 0.69 mK, A' = 3.3 ^iK, C = 50 
AiK, C" = -0.79 /iK, = 2.00 and = 1.93. and 
iJr are, respectively, the longitudinal and transverse compo- 
nents of the magnetic field, and (f) measures the angle between 
Ht and the x axis. Much of the system's dynamics can be 
understood in terms of the double-well potential shown in the 
inset of Fig. |2^, which shows the system's energy as a func- 
tion of the angle between the spin vector and the easy axis (z) 
direction. The spin has 25* + 1 energy levels, which are ap- 
proximate eigenstates of Sz- A longitudinal field tilts the 
potential and at certain values, levels in opposite wells align, 
allowing resonant tunneling between wells that results in a 
marked increase in the rate at which spins flip direction. The 
effect of the fourth and fifth terms in Eq. |2] as well as Ht is 
to break the commutation of T-L and Sz, thereby inducing tun- 
neling. The tunnel splitting between nearly degenerate states 
is readily calculated by diagonalizing Eq.|2] The solid lines in 
Fig. I2J5, for = (mod 90°), show that destructive interfer- 
ence between tunneling paths induces a dramatic suppression 
("quenching") of tunneling at discrete values of Ht for each 
pair of levels ||30|] . (Here and below, we label each state by a 
value TO, its value of {Sz) in the absence of tunneling.) The 
interference effect is largely destroyed when </) is increased 
towards 45° (dashed lines) since Ht then favors one tunnel- 
ing path more than others. The tunneling quenching affects 
which levels are involved in the magnetic relaxation process, 
as evidenced by our data. 

Crystals of Mni2-tBuAc were synthesized according to 
published procedures |27]. A sample was mounted adjacent 
to a Hall sensor that was in turn mounted on a rotator probe. A 
reference sensor on the same chip was used to measure back- 
ground signals. The signals from the two sensors were sub- 
tracted with an analog circuit. We performed extensive mea- 
surements on two samples (A & B). Measurements of sample 
B were performed using a modified apparatus in which the ori- 
entation of both the sample's easy and hard axes relative to the 
field could be adjusted (the latter ex situ). Measurements were 
performed as follows. The sample was rotated so that its easy 
axis was aligned with the external magnetic field, magnetiz- 
ing the sample, i.e. populating one of the wells in Fig.|2]inset. 
Next, the sample was rotated to an orientation that produced 
the desired values of Hz and Ht- The subsequent relaxation 
of the magnetization was monitored as a function of time and 
fit to an exponential decay to extract the associated relaxation 
rate, F. Proceeding this way, we obtained values for the relax- 
ation rate as a function of Hz and Ht and of temperature T 
for relaxation near a tunneling resonance. 

Some data from sample A is shown in Fig. [T] where F is 
plotted as a function of Hz for several values of Ht at T = 
3. 10 K. For each value of Ht, F exhibits a roughly Lorentzian 
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FIG. 1 . (Color online) Measured magnetic relaxation rate as a func- 
tion of longitudinal field near the zero-field resonance for several val- 
ues of transverse field (Ht) from 1.0 kOe to 6.8 kOe in 200 Oe in- 
crements. The data was taken from sample A at 3.10 K. The green 
boxes indicate regions where curves for different values of Ht tend 
to bunch together. The inset shows the classical energy landscape in 
a spherical polar plot for a spin described by Eq[T] The z axis is the 
easy axis (energy minima) while the x and y axes are the hard axes 
(maxima). The value of C has been greatly exaggerated to make the 
four-fold symmetry evident. 



dependence on Hz, peaked at Hz = 0, where tunneling is 
maximum. F generally increases with increasing Ht as the 
tunneling rate is enhanced and, in tandem, the effective en- 
ergy barrier is reduced 113 ill . We note that for some regions 
of Ht and Hz (green boxes) the data are bunched - the relax- 
ation rate changes very little with increasing Ht- The effect is 
much more pronounced in the shoulders of the peak than near 
its center Fig.|2]shows F as a function of Ht for Hz = (up- 
per four data sets, from peak center in Fig.lT]! and Hz = —400 
Oe (lower four data sets, from peak shoulders). All sets show 
a roughly exponential increase in F with Ht- In addition, 
F also exhibits steps and plateaus (the latter corresponding 
to the bunching in Fig. [T]i. These are far more apparent in 
Hz — —400 Oe data. Each step corresponds to a transition 
from one dominant pair of tunneling levels (e.g. m — ±3) 
to another (e.g. to = ±4). Interestingly, theoretical calcula- 
tions of these transitions in which only Ht induces tunneling 
(i.e. when C = C" = in Eq. |2]l predict that the transitions 
should be independent of (/> and much more pronounced on 
resonance (H^, = 0) than when Hz is tuned away from res- 



onance rt32l 13311 . in contrast to our results. Inhomogeneous 



dipole fields can easily wash out any structure near Hz = 0, 
where F depends strongly on Hz - More importantly, the struc- 
ture observed near Hz = —400 Oe (where F is less sensitive 
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FIG. 2. (Color online) (a) Relaxation rate as a function of transverse 
field at Hz ~ and Hz = —400 Oe, taken from the data shown 
in Fig. [T] Some values were obtained via interpolation of the Hz 
dependence. The solid curves are the results of simulations of the 
relaxation rate using Eq.|2]with (j) = 7°, as described in the text, (b) 
Calculated tunnel splitting for energy-level pairs m = ±2, m = ±3, 
and m = ±4, as indicated, for = (solid) and (f} — 45° (dashed). 
The dashed horizontal line running through these calculations marks 
a threshold splitting at which the system makes a transition from one 
dominant pair of tunneling levels to another. The dashed vertical 
lines show the transverse field at which this transition occurs and 
its correspondence to rapid increases in the relaxation rate. Inset: 
Probability current diagram for high-lying energy levels at Ht = 4.8 
kOe. The numerical labels indicate the expectation values of Sz for 
the corresponding energy level. 



to variations in dipole fields) is much more pronounced than 
predicted by the simple model. Including all of the terms in 
Eq.|2]and, especially, the experimentally determined values of 
C and C" induces the tunneling suppression shown in Fig.|2j3. 
These tunnel quenches, in turn, give rise to the steps in F. 

To illustrate this relationship, Fig.|2j3 contains a horizontal 
dotted line - an empirically determined "tunnel threshold". 
When the tunnel splitting for a particular pair of levels ap- 
proaches this threshold, tunneling for that pair begins to be- 
come the dominant relaxation mechanism i3lll32ll34ll35ll . For 
example, at Ht ^ 4.2 kOe (marked by the red vertical dotted 
line), the tunnel splitting for the m — ±3 states reaches the 
threshold and tunneling between these levels begins to dom- 
inate over relaxation through higher lying levels. This addi- 
tional relaxation mechanism produces the rapid increase in F 
near this field. Similar transitions occur when other pairs of 
levels reach the threshold at the values of Ht marked by the 
other dotted vertical lines in the figure. The tunneling sup- 
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FIG. 3. (Color online) Relaxation rate from sample B as a function of 
transverse field for several different values of 0' for Hz = —500 Oe. 
The inset shows the raw relaxation rate while the main figure shows 
the same data after dividing by an exponential function to enhance 
clarity. 



pression effect plays a crucial role here in that it determines 
how rapidly the threshold is crossed as Ht increases. The 
slope of the tunnel splittings curves in Fig. |2]3 is rather steep 
in the vicinity of the threshold right after a quench for = 0. 
In contrast, when = 45°, the tunnel splitting quenches for 
the relevant tunneling pairs are absent and the tunnel splittings 
cross the threshold more gradually, making each transition so 
broad that it overlaps with others and washing out its observ- 
ability. 

This effect is demonstrated in Fig. [3] The inset of the fig- 
ure shows r as a function of Ht for sample B in the vicinity 
of one of the transitions for several values of 0' — + '/'o? 
where 0' is the experimental azimuthal angle of the rotator's 
second stage and 0o is a constant offset representing the ori- 
entation of sample's hard axis. The data features for sample B 
were less distinct than for sample A, possibly because B was 
measured more than a year after synthesis. To enhance clarity, 
we divided the data by gO-Sffr/kOe^ where the coefficient 0.3 
was chosen empirically. The resulting data are presented in 
the main figure. For 0' = 2° and 20°, the sharp transition at 
^ 8 kOe is apparent while it is clearly suppressed for values 
of 0' outside this range. These results suggest 0o ^ 9°, im- 
plying that the hard axes are roughly parallel with the a and b 
crystallographic axes of the rectangulopiped-shaped crystals. 

We also studied the relaxation rate near the n = 1 (H^ ~ 
4.5 kOe) resonance. Data was taken for sample B with 
0' = 2° at values of i/^ ^ 1 kOe below the resonance peak. 
As with the n = {H^ ~ 0) resonance, we see steps in the 
relaxation rate as a function of Ht, as shown in Fig.|4^. (Like 
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FIG. 4. (Color online) (a) Relaxation rate (after dividing by an expo- 
nential to enhance presentation) as a function of Ht for Hz near the 
n = 1 resonance. Data was taken from sample B with experimental 
azimuthal angle (j>' = 2° and set to -1 kOe less than the peak 
of the resonance, (b) Simulations of relaxation rate as a function of 
transverse field using the Giant-Spin Hamiltonian model for SMMs 
with <j) — 7° . The inset shows the double- well potential for the n — 1 
resonance, (c) Tunnel splitting calculations with = 0° for the m = 
-2,3 and m = -3,4 energy-level pairs, as indicated. The dashed hori- 
zontal line represents a tunneling threshold while the dashed vertical 
lines correspond to the transverse fields at which abrupt transitions 
in relaxation rate occur. 



in Fig. [3] we have removed an exponential background for 
clarity.) The steps here also correspond to the transitions be- 
tween dominant tunneling level pairs, as illustrated in Fig.|4j;, 
each occurring when the tunnel splitting for a particular pair 
rapidly crosses a threshold value (horizontal dotted line) in the 
wake of a quench. 

To quantitatively analyze our data, we performed numeri- 
cal calculations of the relaxation rate using a master equation 
approach 1321 l33l 135143711 to treat spin-phonon interactions. 
For calculational ease, we used the spin's energy eigenbasis, 
which incorporates tunneling effects automatically since the 
eigenstates of Eq.|2]are superpositions of Sz eigenstates. We 
neglect off-diagonal elements in the density matrix, a good 
approximation since our experiments were done away from 
the exact resonance conditions where such elements are ap- 
preciable. The master equation governing the population of 
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The phonon transition rates are given by Il32ll38l] 
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where s 



(1) _ 



Ei — Ej, with Ei the energy of level iV (A) — 
f^^A/ksT _ jg j-jjg phonon thermal distribution function, 
p — 1.356 X 10"^ kg/m'^ is the mass density Izsll . and Cs is 
the transverse speed of sound, k'^^ = 1 and k*^^' is a con- 
stant of order unity representing the strength of the associated 
spin-phonon coupling mechanism ||39r. We neglect possible 
collective spin-phonon interactions [4(7]. 

We calculated the relaxation rate by finding the slowest 
non-zero eigenvalue of the rate matrix implicit in Eq. [3] The 



calculated rates are fit to the data in Fig.|2] allowing i 



C and C to be unconstrained parameters. The remaining 
Hamiltonian parameters were fixed and we set — —400 
Oe and (j> — 7°. The results of the fitting are shown by the 
solid curves in Fig.|2^. The calculated rates reproduce the data 
quite well. From the fit, we deduce = 1 122 m/s and k'^^^ = 
1.21. These parameters set the overall scale of the rate and the 
general slope of the rate versus Ht, respectively. They do not 
influence the positions of the steps, which are determined by 
Hamiltonian parameters C and C. The fit yields C = 55 fiK 
and C — —0.81 pK, in good agreement with the values de- 
termined spectroscopically [25]. Using the same parameters, 
we can also calculate the relaxation rate for the n = 1 res- 
onance, shown in Fig. |4j3. Again, the calculations accurately 
reproduce the structure of the measured relaxation rates. We 
do not fit this data because of the reduced quality of the data 
for sample B. 

Our calculations allow us to precisely determine which 
pairs of levels dominate the tunneling process as a function 
of Ht- The inset of Fig|2]3 shows an example of the prob- 
ability "currents" (dashed aiTows) jssll between some of the 
relevant states at one value of Ht {— 4.8 kOe). The state la- 
bels indicate the values of {Sz}- Diagonal arrows correspond 
to tunneling transitions in the energy eigenbasis - such transi- 
tions would be forbidden in the absence of tunneling. For the 
example shown, the states with m = ±3 are clearly the domi- 
nant tunneling levels for this value of Ht- These calculations 
confirm the interpretation of the steps in the relaxation rate 
given above, e.g. for rt = at Ht ^ 4.2 kOe the dominant 
tunneling pair switches from m = ±2 to m = ±3. 

In conclusion, our measurements of the relaxation rate for 
the Mni2-tBuAc SMM as a function of Ht provide the first 
evidence for a geometric-phase interference effect in a truly 
four-fold symmetric SMM. The results also provide clear evi- 
dence that this effect exists in the thermally assisted tunneling 
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regime and allow a clear identification of which levels domi- 
nate the tunneling process. 
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